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Abstract 

We present a measurement of the D°-D mixing parameter yep using a flavor- untagged sample 
of D° — > K^K + K~ decays. The measurement is based on a 673 fb^ 1 data sample recorded by the 
Belle detector at the KEKB asymmetric-energy e + e~ collider. We find yep = (0.21 ± 0.63(stat.) ± 
0.78(syst.) ± 0.01(model))%. 

PACS numbers: 13.25.Ft, 11.30.Er, 12.15.Ff 
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Particle-antiparticle mixing has been observed in several systems of neutral mesons: neu- 
tral kaons, and B s mesons. As in the kaon and 5-meson systems, the D° - D are 
produced in flavor eigenstates. The mixing occurs through weak interactions between the 
quarks and gives rise to two different mass eigenstates |-Di,2 >= p|-D°> ± g|-D°>, where 
p and q are complex coefficients satisfying \p\ 2 + \q\ 2 = 1. The time evolution of flavor 

eigenstates, D° and D° , is governed by the mixing parameters x = (mi — 7712) /Y and 
y = (Yi — Y 2 )/2Y, where and are the masses and widths of the two mass eigen- 
states D12, and T = (Yi + Y 2 )/2. In the Standard Model (SM), D° - D mixing is strongly 
GIM suppressed for d and s quarks and CKM suppressed for b quark box diagrams, and is 
dominated by long distance effects As the mixing rate is expected to be small within 
the SM, it is sensitive to the contribution of new, as yet unobserved processes and particles. 
The largest SM predictions for the parameters x and y, which include the impact of long 
distance dynamics, are of order 1% [l|. Various D° decay modes have been used to measure 
or constrain x and y 0]. Evidence for D°-D^ has been found in D° — > K + K~ /7i + 7i~ 0, 0], 
D° -> K + 7T~ @, Q and D° -> K+tt'tt Q decays. The world average @, Q of D° mixing 
parameter yep measured in D° — > K + K~ /it + it~ decays is yep = (1.132 ± 0.266)%, where 
Ucp — V if CP is conserved. Here we study the self-conjugate decay D° — > K^K + K~ [10]. 
The time dependent decay rate of an initially produced D° or D° can be expressed as 

HI M 
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T 

where r = 1/Y is the D° lifetime, sq and s + are invariant masses squared of K + K~ and 
K$K + pairs, respectively. The decay amplitudes A\ and A 2 can be expressed with D° and 
D° decay amplitudes A and A as 

A(s ,s+) = 5> r .e^A.(so,s+) (3) 

r 

A(s ,s + ) = ^a r e^A(s ,s+) (4) 

r 

A\(so, s + ) = - (A(s , s + ) + A(so, = ^2 CP = +1 and flavor eigenstates (5) 
A 2 (so, s +) = 2 ("4( s o, s +) — A(s , s + )^j = CP = —1 and flavor eigenstates, (6) 

where A and A are summed over resonant contributions r found in D° — > K®K + K~ de- 
cays. In the limit of CP conservation a r = a^, 4> r = (fi r and A(s , s + ) = A(s ,s_). The 
existing Dalitz plot analyses of D° — > i^f^+f^decays [l3|, [13] observed contribution of 
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FIG. 1: so (left) and s + (right) Dalitz plot projections of \M(s ,s + )\ 2 (black line), 
|-4i(so> s +)| 2 (blue line ), (^(so, s + )\ 2 (red line), 2Re\A-\ (so, s+)A%(sp, s+)\] (green dotted line) 
and 2Im[Ai(so, s + )A2(sq, s+)\] (green dashed line) for Dalitz model given in [14]. 



CP (iv%(980)°, iv°0(lO2O), ^/ (1370), AT°/ 2 (1270), AT°a (1450)°, K° s f Q (980)), Cabbibo- 
allowed (K~a (980) + , K~a (1450) + ) and doubly Cabbibo-suppressed (K+a (980)") fla- 
vor eigenstates. Figure CD shows time integrated sq and s + projections of |.M(so, s + )\ 2 , 
\Ai(s ,s+)\ 2 , \A 2 {s ,s + )\ 2 , 2Re[A 1 {s ,s+)A2{s ,s + )\] and 2Im[Ai{s , s + )A* 2 {s , s+)\] ob- 
tained by Dalitz model given in Ref. |l4j . The integral of 2Re[A\(s , s + )A%(s , s + )\] and 
2Im[Ai(so, s + )Al(so, s+)\] over s + yields 0. 

The |^4i| 2 and l^] 2 parts of the decay rate have different time dependence (Eq. [1] and 
[2]) and also very different dependence in the sp (Fig. [1] (left)). In any given sp region the 
lifetime of D° candidates is given by 

^' = h— — + 0-fi) T ^ ! —, (7) 

i + yep i - yep 

where r is the mean D° lifetime 1/r, f\ = §\Ai\ 2 / /(|Al| 2 + l^l 2 ) and CP conserva- 
tion is assumed. The lifetime difference of D° candidates in two different regions is then 
proportional to the mixing parameter yep 

*r = ^ = *»> 1 + y J^J)„ _ fl) « VcAf! ~ /!). (8) 

The best m(K + K~~) intervals from which D° lifetimes are measured and compared are 
those that minimize the statistical uncertainty on yep and are found to be: region 
around 0(1020) peak m(K + K~) e [1.015, 1.025] GeV/c 2 (denoted as ON) and intervals 
m{K+K-) G [2m K ±, 1.010] GeV/c 2 and m(K + K~) e [1.033,1.100] GeV/c 2 (the union of 
this two intervals is denoted as OFF), where m^± is the nominal K mass. 

The data were recorded by the Belle detector at the KEKB asymmetric-energy collider 



151 ] . The Belle detector is a large-solid-angle magnetic spectrometer that consists of a silicon 
vertex detector (SVD), a 50-layer central drift chamber (CDC), an array of aerogel threshold 
Cherenkov counters (ACC), a barrel-like arrangement of time-of-flight scintillation counters 
(TOF), and an electromagnetic calorimeter (ECL) comprised of CsI(Tl) crystals located 
inside a superconducting solenoid coil that provides a 1.5 T magnetic field. An iron flux- 
return located outside of the coil is instrumented to detect K\ mesons and to identify muons 
(KLM). The detector is described in detail elsewhere [l6|. Two inner detector configurations 
were used. A 2.0 cm beampipe and a 3-layer silicon vertex detector was used for the first 
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FIG. 2: The distribution of m(K s ) with m{K s K+K~) £ [1.85,1.88] GeV/c 2 (left) and 
m(K s K + K~) with m(Ks) G [0.490,0.505] GeV/c 2 (right). Superimposed on the data (points 
with error bars) are projections of the m(Ks) — m(K s K + K~) fit (result from the fit (solid blue 
lines) , signal contribution (solid green line) , true K s (solid black line) and rest of the background 
(solid red line)). 

sample of 156 fb -1 , while a 1.5 cm beampipe, a 4-layer silicon detector and a small-cell inner 
drift chamber were used to record the remaining 517 fb _1 of data. 

The candidates are reconstructed in the 7r + 7r~ final state; we require that the pion 
candidates form a common vertex at least 0.9 mm from the e + e~ interaction point (IP) in 
plane perpendicular to the beam axis and have an invariant mass within ±30 MeV/c 2 of 
K s nominal mass. We reconstruct D° candidates by combining the K s candidate with two 
oppositely charged tracks assigned as kaons. These tracks are required to have at least one 
SVD hit in both r — cf) and z coordinates. A D° momentum greater than 2.55 GeV/c in 
the e + e~ center-of-mass (CM) frame is required to reject D mesons produced in B mesons 
decays and to suppress combinatorial background. 

The decay point of D° candidate is determined by refitting one of the charged kaons and 
K s candidate to a common vertex [l8|; confidence levels exceeding 10~ 3 are required for 
the both fits. Out of two possibilities the one with lowest \ 2 value of the fit is used. In 
addition we require that K S K + K~ and K + K~ combinations originate from the common 
vertex by rejecting candidates of this two fits with confidence levels lower than 10~ 3 . The 
D° production point is taken to be the intersection of the D° momentum vector with the 
IP. The proper decay time of the D° candidate is then calculated from the projection of 
the vector joining the production and decay points, L, onto the D° momentum vector, 
t = (mpo /po)L ■ {pd/pd), where m^o is the nominal D° mass. The decay time uncertainty 
a t is evaluated event-by-event, and we require at < 600 fs (the maximum of o t distribution 
is at ~ 230 fs). 

The signal and background yields are determined from a two-dimensional fit to the in- 
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variant masses of K s and D° candidates. According to Monte Carlo (MC) simulated distri- 
butions of m(Kg) and m(K s K + K~), events can be divided into three categories: (1) signal 
D° — > K S K + K~ decays; (2) true K s candidates combined with random charged kaons (one 
or both); and (3) rest of the background. We parametrize the signal shape by a sum of 
a three two-dimensional-Gaussian function and a product of two one-dimensional-Gaussian 
functions used to describe long tails in both variables (the contribution of the latter is small, 
~ 0.1%). The second category is described by a sum of three Gaussians for m(Kg) and a 
linear function for m(K s K + K~). The third category is described by a product of the linear 
functions. The m(K s ) and m(K s K + K~) distributions are shown in Fig. [2] along with the 
projections of the fit result. The MC simulation shows that a small fraction (~ 0.1%) of 
events are decays of D° mesons to K + K~tt + tt~ final state (charged pions do not originate 
from Kg decay). These events are peaking in m(K s K + K~), but not in m(K s ). The projec- 
tions of m(K s K + K~) for events in m(K s ) sidebands are checked for possible contribution 
of D° —>■ K + K~7t + 7t~ decays. We find no contribution of D° — > K + K~7t + tt~ decays. The 
fit is performed to obtain scaling factors for the background fractions, and then tune them 
in the MC event-by-event in order to achieve better agreement in m(K s ) and m(K s K + K~) 
distributions between MC and data events. 

The sample of events for the lifetime measurement is selected using \m'(Kg)\ and 
\m'(K s K + K~)\, where m'(K s ) and m! (K S K + K~~) are rotated K s and D° candidate masses 
according to 

where m K o = 497.57±0.01 MeV/c 2 and m D o = 1864.96 ±0.01 MeV/c 2 are fitted Kg and D° 
masses, a(K° s ) = 1.826±0.006 MeV/c 2 and a(K s K + K~) = 2.915±0.009 MeV/c 2 are widths 
of the core Gaussian function and p = 0.602 ±0.002 is the correlation coefficient. The above 
uncertainties are statistical only. We define the signal box in the plane of rotated masses 
m'(Kg) and m'(K s K + K~) in order to minimize correlations. Signal window in \m'(Kg)\ 
and \m'(K s K + K~)\ is chosen to minimize the expected statistical error on yep, using the 
tuned MC: we require \m'(Kg)\ < 3.9 and \m'(K s K + K~)\ < 2.2. The selection criteria on 
a t and Kg candidate flight distance in r — <fi plane, given above, are determined in the same 
way. We find 139 x 10 3 signal events with purity of 94%. 

The lifetime difference A r (Eq. [8]) is determined from D° — > K S K + K~ proper decay 
time distributions by measuring lifetime of signal events in ON and OFF m(K + K~) regions. 
The lifetime of signal events is obtained in the following way. For each event category i the 
proper decay time distribution Pi(t) is assumed to be either exponential or a delta function, 
convoluted with a resolution function Ri(t). The distribution for all event categories is then 

P(t) = J Ep i P i (t)®R i (t), (11) 

i 

where Pi — iVj/ Y^j Nj is a fraction of the category i. By grouping the events into the signal 
and background one can also write 

P{t) = p-e- 1 ^ <g> R s (t) ± (1 - p)B(t), (12) 
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m{K+K-) 


^sw N sh <i> sw [fs] <t> sb [fs] p[%] 


T s + 1 [fs] 


OFF left 

ON 
OFF right 


19618 763 400.2 ±4.5 121.2 ± 27.7 96.11 ± 0.14 
66112 2104 403.0 ± 2.4 41.2 ± 13.8 96.82 ± 0.07 
40634 4879 381.6 ± 3.2 138.6 ± 10.2 87.99 ± 0.16 


411.5 ±4.8 
414.9 ±2.6 
414.7 ±3.9 



TABLE I: Numbers of events in the signal window N sw and sideband iV s b, mean proper decay 
times of events in the signal window < t > sw and < t > s b, fraction of signal events in the signal 
window p = 1 — iV s b/iV sw and reconstructed lifetime r s + to (Eq. [T4"j) shifted for resolution function 
offset obtained on untagged real data sample. 



where the first term represents the measured distribution of a signal with lifetime t s , R s (t) 
is a signal resolution function and p = N s / (N s + ]V&) is a fraction of signal events. The last 
term represents the distribution of background events. The mean of the above distribution 
(Eq. [T2D is 

<t>=p(T 8 + t ) + (l-p) <t> b , (13) 

where to is the mean of the signal resolution function R s (t) and <t> b is the mean lifetime 
of the background. The lifetime of signal events, shifted for the resolution function offset, 
can be calculated from Eq. 



with uncertainty 



r <t> -(1 -p) <t> b , , 

T s + t = (14) 

p 



(15) 




where a, a b and a p are determined from the proper decay time distributions of all events 
P(t) and background events B(t) in the following way 



rms(P) rms(B) p(l—p) 

a = -=— , a b = — -=- and a p - ■' 



N VN b y N 

The B(t) distribution of background events populating the signal window is approximated 
by the proper decay time distribution of events taken from m'(KgK + K~) sideband of equal 
size as signal window. No scaling factor is needed, since the background events are linearly 
distributed in m'(KgK + K~). The tuned MC is used to select the sideband region that best 
reproduces the timing distribution of background events in m'(KgK + K~) signal window, 
which is chosen to be 9.7 < \m' \K Q S K + K')\ < 11.9. 

In Table [U the numbers of reconstructed events in the signal window iV sw and sideband 
N sh , mean proper decay times of events in the signal window <t> sw and <t> sb , fraction of 
signal events in the signal window p = 1 — N sh /N sw and reconstructed lifetime t s + t (Eq. 
[T4l shifted for resolution function offset obtained on real data sample are given for 3 different 
regions: OFF left {m{K + K~) < 1.010 GeV/c 2 ), ON (1.015 < m{K + K') < 1.025 GeV/c 2 ) 
and OFF right (1.033 < m(K + K~) < 1.100 GeV/c 2 ). Figure [3] shows proper decay time 
distributions for events populating OFF left, ON and OFF right m(K + K~) regions. 
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FIG. 3: Proper decay time distributions for events populating m(K + K ) < 1.010 GeV/c 2 (left), 
1.015 < m(K + K~) < 1.025 GeV/c 2 (middle) and 1.033 < m(K + K~) < 1.100 GeV/c 2 (right). The 
hatched area histograms show the contribution of events populating the m' {K^K + K~) sideband. 

To obtain yep from measured A r (Eq. EJ) the fraction difference /° N — / 1 OFF is needed. 



Dalitz models of D° — > K$K + K decays given in Ref. [13, ll4j are used to fit the Sq 



distribution. The s distribution of signal events is parametrized as 

P sig = e(s ) J e(s+)\M(s , s + )\ 2 ds + , (16) 

where \J\A(s , s + )\ 2 is the time integrated decay rate (Eq. [I]), and e(sq) (e(s + )) is the 
reconstruction efficiency in so (s+) determined from a sample of MC events in which the 
decay mode was generated according to phase space. Efficiency in so and s + is assumed 
to be factorizable. No significant difference is observed between e(s + ) obtained for events 
populating ON and OFF sq regions. All phases, amplitudes, masses and widths of the 



resonances are fixed to the values determined in Ref. [131 . [14J , except for the amplitudes of 
#§0(1020) and #-a (980)+ (#-a (1450)+) contributions using model from Ref. (Q). 
The free parameters of the fit are also the coupling constant gxx of a coupled channel BW 
[Tol which describes the ao(980) resonance and the mass and width of the 0(1020) resonance 
in order to account for mass resolution effects. To describe background events in the Sq 
distribution, events from the m'(KgK + K~) sideband are taken. The \ 2 test of the MC 
So distributions of background events taken from the signal window and sideband yields 
X 2 /ndf = 136/99. The fraction of signal events in the signal window p = 1 — N 8 \ ) /N sw is 
determined from the numbers of events in the signal window iV sw and sideband iV s b and it is 
fixed parameter of the fit. Figure H] shows the fit result to the so distribution for the Dalitz 
model given in Ref. [13]. The \ 2 l n df value of the fit is 431.8/230 using the Dalitz model 
from Ref. [l3[ and 291.7/230 using the Dalitz model from Ref. [13]. In Table HT1 fractions 
/° N and fi F¥ and the fraction difference f± N — / 1 OFF are given for both Dalitz models. 
Although the models are very different, with different resonant structure lj|, the fraction 
differences calculated for each model are in agreement. 

The reconstructed lifetimes shifted for the resolution function offset, r s + t , of D° candi- 
dates in ON and OFF regions are 414.9 ± 2.6 fs and 413.6 ±3.1 fs, respectively, from which 
A T = (—0.16 ± 0.48)% is obtained. We assumed that the resolution function offset, to, is 
equal for the events populating the ON and OFF regions and much smaller than D° lifetime. 
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FIG. 4: s distribution of D° 
given in Ref. 



K decays with superimposed fit results with Dalitz model 



14i | (right). The blue solid line is the overall fitted function and the red line is the 



background contribution. 





Nominal 


Fitted 


Model 


jON jOFF ^ON _ yOFF 


fON f OFF fON f OFF 
Jl Jl Jl Jl 


4 res. [13] 
8 res. [14] 


0.117 0.847 -0.730 ±0.031 
0.124 0.877 -0.753 ±0.004 


0.113 0.844 -0.732 ±0.003 
0.111 0.880 -0.769 ± 0.005 



TABLE II: Fractions ff" and /f FF (f° N l° FF = ^qn /off \M 
fraction difference /f N — /P FF for the two Dalitz models Ref. 



/ 



'On/off(I--4i| 2 + l^2| 2 ) and the 
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141 ] . The nominal values are 



calculated using the given Dalitz models in Ref. [13, [Uj and fitted values using the obtained values 
of free parameters of the fit to the sq distribution. Uncertainties on /f N — /P FF were calculated 
using the statistical errors of amplitudes and phases given for each model, without taking into 
account any correlation between the amplitudes and phases. 



Using the Eq. [8] and the fraction difference f\ 



ON 



distribution with Dalitz model given in Ref. [14 



_ yOFF _ 

yields y C p 



-0.769, obtained by fitting s 
-- (0.21 ±0.63(stat.) 



o. 



We consider systematic uncertainties arising from both experimental sources and from 
the D° — > KgK + K~ model. First, we check on the MC sample if the resolution function 
offsets, £q and t® FF are equal. They are in agreement within the statistical uncertainty 
and small (t = 0.7% ■ t d o). Next, we vary the sideband in m'(KgK + K~) used to describe 
the background populating the signal window and measure for each sideband the A T . For 
different sidebands used the obtained A r values are in agreement. The maximal difference 
in A r was taken to estimate the systematic uncertainty. Finally, possible systematic effects 
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SvKtpiTiaflP PTTOr ( % 1 


Resolution function offset difference t[j )FF — t$ 


±0.21 


Selection of m'(KgK + K~) sideband 


±0.35 


Variation of selection criteria 


±0.44 


Total 


±0.60 



TABLE III: Sources of the systematic uncertainty for A T . 



of selection criteria were studied by varying the signal box sizes, and cut values on a t and 
Kg flight distance in r — <p plane. Again no statistical significant deviation was observed 
and the maximal difference in A r was taken to estimate the systematic uncertainty. We add 
all different sources in quadrature to obtain the overall experimental systematic uncertainty 
summarized in Table [TTT1 

The systematic uncertainty due to our choice of D° — > K$K + K~ decay model is eval- 
uated as follows. First, we compare the fraction difference — /° FF obtained using the 



Dalitz Models in Ref. [13l. [14j . Despite the differences between the two models in terms 
of the resonant structure 19J, the fraction differences / 1 ON — f® FF (given in Tab. [TT]) are 
in agreement. We assign 3% relative error for measured yep due to small difference in the 
above fractions. An additional 2% relative error for measured yep is assigned due to the 
small difference between fitted and nominal values of fraction difference — / 1 OFF (given in 
Tab. [XT])- The real and imaginary part of the interference term AiAlj in the decay rate (Eq. 
[T]) are zero after integrating over the s + . Since the reconstruction efficiency is not constant 
in s + , this is not entirely true. However, even if the observed s + reconstruction efficiency 
is taken into account this has negligible effect and Eq. [8] still holds. This was also verified 
by MC with non-zero x and y values of mixing parameters, where the detector response 
was simply simulated by randomly rejecting events according to the observed dependence 
of efficiency in s+. The difference between the obtained A T values (with and without tak- 
ing into account the efficiency in s + ) are in agreement within statistical uncertainty, so no 
additional systematical uncertainty is assigned. Adding all variations in quadrature, the 
obtained relative model systematic uncertainty is 4%. 

In summary, we determine yep by measuring the difference in lifetimes between D° 
mesons decaying to K^K + K~ in two different m(K + K~) regions with different contributions 
of CP even and odd eigenstates to be 

y C p = (0.21 ± 0.63(stat.) ± 0.78(syst.) ± 0.01(model))%. 

The result is in agreement with world average of yc p of previous measurements 0, [9| . 
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